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CN ■ We introduce two simple models which feature an Alice electrodynamics phase. In 

a well defined sense the Alice flux solutions we obtain in these models obey first 
■ order equations similar to those of the Nielsen-Olesen fluxtube Q in the abelian 

higgs model in the Bogomol'nyi limit. Some numerical solutions are presented as 
well. 



a: 
> 

X 1 Alice electrodynamics 

a i 

Alice Electrodynamics is a theory with gauge group H = U(l) x Z 2 , i.e. an absolute 
minimally non-abelian extension of ordinary electrodynamics, where charge conjugation 
has been turned into a local symmetry. As this non-abelian extension is discrete, it 
only effects 'electrodynamics' through certain global (topological) features, such as the 
existence of Alice fluxes and Cheshire charges @, j3|. In an Alice phase of some theory one 
has the possibility of a topological stable flux tube or string coexisting with an unbroken 
U(l), because the connectivity of H is nontrivial: U (H) = Z 2 . The Alice phase is 
usually obtained by spontaneous breaking of a larger, continuous non-abelian symmetry 
group. The original Alice model studied by Schwarz 0, is a SU{2) theory spontaneously 
broken down to a U(l) x Z 2 by a higgs field in the 5— dimensional representation of 
(see alsoQ). The higgs field is chosen in this representation, because it is the smallest 
irreducible representation which admits H as a residual symmetry group and allows for a 
single valued vacuum configuration that supports Alice fluxes. 
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In this paper we will discuss two alternative models, which support an Alice phase. Before 
doing so, we briefly review the salient features of the model discussed in ||, [§. 
The action is given by: 

S = Jd 4 x {Tr {--F^Fp, - l -D^D^ } + V(®)} , (1) 

where the higgs field $ = $ ab is a real symmetric traceless 3x3 matrix. 
The most general renormalisable potential is given by: 

V = -^fi 2 Tr$ 2 - ^ 7 Tr$ 3 + ^A(Tr$ 2 ) 2 (2) 

By a suitable choice of parameters the higgs field will acquire a vacuum expectation value, 
$o- In a gauge where $o is diagonal it takes the form $o — diag(a, b, —a — b). For a certain 
range of potential parameters one furthermore has that a = b, so that $o is given by: 



$o = a (3) 




7± \J 7 2 +24/1 2 A 

with a = b = — v 12A . Indeed, this ground state is invariant under rotations around 

the T3-axis (U(l)) and invariant under rotations by an angle it around any axis perpen- 
dicular to the T3-direction (Z2). These two transformations do not commute with each 
other, in fact they anti-commute, so the resulting residual gauge group is indeed U (1) ix Z 2 . 
This means that we have Alice electrodynamics as the low energy effective theory in this 
model. 

An alternative way to see the structure of the residual gauge group, is to think of the 
higgs field as the symmetric traceless product of two vectors, 4>i an d 4>2, 



= + (4) 

If both isovectors, <pi, are non zero, there is in general only a Z 2 gauge symmetry left, 
4>i ~^ 4>'i = ~ ( t ) i- However, in case that both isovectors are (anti-)parallel, the gauge group 
is U(l) ix Z 2 . If one of the isovectors is zero, the gauge group is not broken at all and the 
symmetry remains SU(2). These are the residual gauge groups which one may encounter 
in this model. It is easy to show that the case where the two isovectors are (anti-)parallel, 
corresponds to the situation where $ = $ • 



1.1 The Alice flux solution 

In this section we will present explicit regular solutions, corresponding to an Alice flux 
tube along the z-axis, which where constructed in ||. 

To have a static finite energy solution, all terms in the energy density should go to zero 
at spatial infinity. Thus the covariant derivatives need to vanish at spatial infinity. Let's 
look at the angular derivative, the condition Dg§ = tells us that the higgs field has the 
following form at spatial infinity. 

$(0) = ^(^(O)^- 1 ^) (5) 
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with 



S(9) = exp{e I rA e dO} (6) 
Jo 

Since we are looking for solutions which correspond to an Alice flux, S(27r) needs to be 
an element of the disconnected part of the (residual) gauge group. A simple choice for A$ 
doing this is A e = ^T x . 
This leads to the ansatz: 

$(r, 9) = e~$(r)e ~ (8) 
where the tensor $(r) is conveniently parameterized as, 

(I \ 

$(r)=m(r) -| + q{r) I § I (9) 
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The part proportional to m(r) is the part of the higgs field that is invariant under rotations 
generated by T\. The boundary condition at spatial infinity is m(oo) = q(oo), implying 
that $(oo) is of the form (||), i.e. the residual symmetry is U(l) x Z2 indeed, where the 
electrodynamic U(l) is generated by T3. At the origin, m and q have to satisfy different 
boundary conditions; the field q(r) needs to go to zero. The term proportional to m(r) 
is invariant under 7\ rotations, therefore m(r) does not need to go to zero. Again, this 
means that the higgs field is of the form (|3p, i.e. the unbroken gauge group is U(l) x Z 2 . 
However, the unbroken U(l) is generated by T\. Finally, the field a(r) needs to be zero 
at the origin and unity at spatial infinity. 

Inserting this ansatz in the field equations gives, after suitable rescalings, the following 
set of equations. 

d 2 r a{r) - -d r a(r) = 9g 2 (r)(a(r) - 1) (10) 



r 



i ;l , _ • ' -_QV' 



df,q(r) + -d r q(r) = K w ' y ' + £(9q 2 (r) + 3m 2 (r) - 2)q(r) 



+2xm(r)q(r) (11) 
d 2 r .m(r) + -d r m(r) = £(9g 2 (r) + 3m 2 (r) - 2)m(r) + x(3g 2 (r) - m 2 (r)) (12) 

We summarize the boundary values for the rescaled fields below: 



field r —y r^ 00 

a(r) 1 

q(r) 9(00) 

m(r) constant m(oo) 



where, 



- X ± y/ X 2 + gg _ 

m(oo) = g(oo) = — = a(£, x) (13) 
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with £ = -4 and x 
The system 



fie 2 



12|) was solved numerically with the help of a relaxation method in g. 
The solution for the potential parameter values £ = 1 and x = — 1 is given in figure [1] . 
The situation at hand is reminiscent to the one considered by Witten || for a U(l) X £7(1) 
model, in the sense that we have an unbroken [/(l) XZ2 at the core and a different U(l) KZ2 
at infinity. However, the crucial difference is, that our 'two' Z7 (1) gauge groups do not 
commute with each other. 




Figure 1: A regular solution for the fields with an Alice flux for £ = 1 and x = — 1. 



Interestingly, there is also another solution to the field equations, which we briefly like to 
discuss. If x = there is a solution with m(r) = 0. After a rescaling of q(r) one finds 
exactly the same equations as were obtained in the Nielsen Olesen (NO) model by [[F] for 
the minimal flux n = 1, provided we set the value of A = 2£. Numerical solutions to these 
equations have been studied before. For a special value of A one obtains the solutions 
by solving, Bogomol'nyi type, first order equations, signaling the possibility of extending 
the model to a super-symmetric one whose super-symmetry increases for this value of the 
parameter. The residual symmetry of this solution in our model is Z 2 . One may wonder 
whether in our model this is a stable solution. In the case of 7 = 0, i.e. x = 0> the 
potential (0) has the form: 

V = -fi 2 X + XX 2 (14) 



The minimum of this potential is obviously given by X = |r, with X 
in the components m and q this gives: 



cim 2 + c 2 q 2 
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lTr($) 2 . Written 



(15) 



with Ci, C2 > 0. A simple rescaling of m and q yields: m + q = 1. As we require finite 
energy, this is one of the boundary conditions for the fields m and q at spatial infinity. 
We now see that the boundary condition of the new solution, m = 0, may be continuously 
changed to one where q = 0. But then the higgs field does no longer stabilize the flux, 
which means that the flux will decay by spreading out and losing more and more energy. 
Through this process we end up in a "rotated" Alice phase of the theory where the U(l) 
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generator points in the internal direction of the flux we started of with. When going from 
the m = to the q = boundary condition we pass an Alice phase whose £7(1) generator 
is perpendicular to the internal direction of the flux we started with. The upshot of these 
observations is that, if one wants to have a stable Alice flux in an Alice phase, one needs 
to have 7 / 0. Or, stated the other way around, if we do have 7 = and are in an 
Alice phase with an (unstable) Alice flux, the Alice flux will decay and change the Alice 
phase into another Alice phase, whose £7(1) generator is in the internal direction of the 
flux we started of with. This concludes what we have to say about conventional Alice 
electrodynamics, in the remaining sections of the paper we will focus on some alternative 
Alice models. 



2 Alternative Alice models. 

In this section we introduce two alternative models, which exhibit an Alice electrody- 
namic phase. In these alternative models we choose the higgs field(s) in the adjoint (3) 
representation of 50(3). This obviously means that the higgs field is not single valued in 
the presence of an Alice flux, but this can be "solved" in two more or less similar ways. 
One way is to put the internal space of the higgs field (A) on a Z 2 orbifold, i.e. X and 
— X are identified with each other. The other way is to use (at least) two higgs fields and 
put the total internal space of these two higgs fields (X and Y) on a S% orbifold, i.e. you 
identify the points (X,Y) and (Y, X). 
The action we use for both models is given by: 

d 4 x {TV {--F^F^ - -D^XD^X - -D^YD^Y + |[X,lf } 

+^{Tr{X 2 + r 2 }-/ 2 } 2 } (16) 

Both theories allow the presence of an Alice flux. In the S2 model it means that one 
studies the twisted sector of the theory. 



2.1 Alice flux solutions 

We now turn to the construction of regular cylindrically symmetric (numerical) solutions 
corresponding to an Alice flux. At spatial infinity one has DgX = 0, implying that the 
higgs field should have the following form at spatial infinity. 

X{0) = S{6)X(0)S- 1 (e) (17) 

with S(9) the same as in @. The flux associated with S(9), is topologically stable if the 
element of the gauge group associated with S(2tv) is an element of the disconnected part 
of the residual gauge group. A simple choice [] is: A e = -^-T\. This puts the Alice flux in 
the internal T\ direction. Writing the higgs field as X = x a T a , it follows that X{9) has 

x At this point we can not yet say that this is an element of the disconnected part of the residual gauge 
group, but this will be done consistently below. 
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the following form: 




X{9)-- 

Thus you see that X(2ir) is given by: 

/ x 1 (0) 
X(2n) = -x 2 {0) 
V -*s(0) 

The same holds, of course, for the other higgs field. Because the two models differ slightly 
in constructing the ansatz, we will treat them separately for the moment . 

The Z 2 model: 

For the Z 2 model the boundary condition specified above implies that either x\ = or 
%2 — x 3 — 0. Only in the first case, however, is S(2iv) an element of the disconnected 
part of the gauge group. Thus we have to put x\ = 0. Later we will see that this choice 
is important in order to obtain first order equations. At this point it is convenient to 
introduce a different basis for the generators of the gauge group, a basis naturally linked 
to the orientation of the higgs field. Its elements are given by: 

S a (8) = e^T a e^~ (18) 

Now we write the higgs field as X = x a S a , where also in this language one has to put 
Xi — to secure the possibility of a topological stable solution. 

In this model a single higgs field would suffice, but for reasons of similarity we will use 
two. Our ansatz than reads: 

(19) 

(20) 
(21) 









2er 


X 


= a(r)S 3 


Y 


= c(r)Si 



The ^2 model: 

The 'double valuedness' is only allowed if one uses an orbifold interpretation. So we im- 
pose a strict relation between X and Y. 

( Xl (0) \ ( yi {9) \ 
X(6 + 2n)= -x 2 (8) = y 2 {6) =Y(6) 
V -x 3 (8) J V 2/3(0) / 

Leading to: x\ = yi , x 2 — —y 2 and x 3 = —yz- Again we are going to work with the 
twisted generators (18). A consistent ansatz is the following one: 

A e = ^-Sx (22) 
2er 

X = a(r)S 3 + c(r)Si (23) 

Y = -a(r)S 3 + c(r)Si (24) 
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For both cases one may insert the appropriate ansatz in the field equations. This yields 
after a suitable rescaling, the same set of differential equations for both models: 



dlair) - -d r a(r) = (a(r) - l)a 2 (r) (25) 
r 

dVr) + -d r a(r) = ( Q ( r ) ~ X ) a(r) + Xa(r)(a 2 (r) + c 2 (r) - 1) 

+7c 2 (r)a(r) (26) 

<9 2 c(r) + -<9 r c(r) = 7a 2 (r)c(r) + Ac(r)(a 2 (r) + c 2 (r) - 1) (27) 
r 

The asymptotic values of the fields are are as follows: 



field 


r -> 


r^ oo 


a{r) 





1 


a(r) 





1 


c(r) 


constant 






The boundary conditions are such that S(27r) is an element of the disconnected part of 
the residual gauge group. 

We have constructed numerical solutions to these equations, for different values of A (and 
7), with the use of a "shooting" method, see figures H and || As a matter of fact, we 
only found solutions for which c(r) = 0, although our starting values were chosen quite 
general. This implies that there is no dependence of the solutions we found, on 7. 

In fact if c(r) = 0, the equations become the same as in the case of a NO flux with 
the critical value of the Landau coupling parameter leading to first order Bogomol'nyi 
equations. However, there is an important difference with the NO case. The "winding" 
number of the Alice flux is fractional and equals value which is not admissible 

in the NO model. This is clearly a consequence of the different breaking schemes of the 
theories in question. 

There is a special role in these theories for the parameter 7, if we set 7 = the equations 
are very similar to the equations ([10|)-(|l2") with x = . Though 7 appears to play no 



role as long as c(r) = 0, this is not quite the case. We don't want 7 to vanish because 
than we run more or less into the same problem as in the conventional model for Alice 
electrodynamics with x = 0. The solution with c(r) = would still be a solution of the 
field equations, but the flux would no longer be stable. It would be allowed to decay into 
the vacuum. In fact in the alternative models it is quite clear what happens at 7 = 0, 
the potential term proportional to 7 - assuming it is nonzero - ensures that there is no 
continuous path in the vacuum manifold connecting the c = to the a = boundary 
condition^. If 7 = such a path does exist. 

There is a simple relation between the Z2 model with one higgs field and the S2 model 
with two higgs fields, in the presence of an Alice string. In the presence of an Alice string 
the field component of the higgs field parallel to the Alice flux is zero in the Z 2 model, 
whereas in the £2 model this is in general only true far away from the core. So, in some 
sense the Z2 model is a long wavelength approximation of the S2 model, but remarkably 



2 Rcmcmber, there is also the boundary condition a 2 + c 2 = 1. 
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enough, it does support solutions which are regular everywhere nevertheless. The action 
of both models becomes the same, up to a rescaling, if the components parallel to the 
Alice flux, of the higgs fields in the S 2 model, are set equal to zero. 




2.2 First order equations 

As mentioned before, if one sets c(r) = 0, the set of equations, ( f25H2T| ), reduces to the 
same set that one would obtain in the NO model for a solution with winding number 
n = \ . It thus appears that one can, in the sector that contains a topologically stable 
Alice flux, project both theories on a sector of the NO model. This raises the question 
whether it would be possible to find first order equations in both models. In the I2 model, 
with only a single higgs field, this projection is the clearest. For the rest of this section 
we will therefore concentrate on this case. 
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One of the features of the NO theory is that for a certain value of the coupling constant 
A, the solutions can be obtained from first order equations. These first order equations 
can be found a la Bogomol'nyi, by rewriting the energy density as a sum of squares plus 
a topological term. In the case of static solutions, the energy density and the Lagrangian 
differ only by a sign, implying that extrema of the energy are also extrema of the La- 
grangian. Consequently, solutions of minimal energy are stable static solutions of the full 
set of second order field equations. The energy of the Z 2 model is given by: 

E = \f d 3 x{ Tr( Ef + B 2 + {D t Xf + (D t X) 2 ) + ^( (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 - 1) } (28) 

In the static case with no electric fields, in the gauge A t = 0, one has Ei = d t Ai = and 
D t X = d t X = 0, reducing the expression to, 

E=\j d 3 x{ Tr( B 2 + (D t X) 2 ) + ^((x 1 ) 2 + (a; 2 ) 2 + (x 3 ) 2 - 1) } (29) 

Restricting ourselves to the plane, the energy written in components is given by: 

E = \ J d 2 x{{B]) 2 + (B 2 ) 2 + {Blf + {{D v Xff + ((ZVO 2 ) 2 + {{D v Xff 

+±((xr + (x 2 ) 2 + (x 3 ) 2 -m (30) 

where now the upper label refers to the internal directions (in the normal, non-twisted 
basis) and the lower label to the spatial directions. From this expression we are unable to 
obtain first order equations, however, if we restrict ourselves to the subspace of solutions 
containing an Alice flux, there is something we can do. Let's call the internal direction 
in which the Alice flux 'points' the T' direction, so if one is looking for topological stable 
fluxes one needs to have X\ = 0, as argued before. In that case we may write the energy 
as, 

E = \ ld 2 x{(Bl) 2 + (Bl) 2 + (B!) 2 + ([A 1 ,X] 1 f + ([A 2 ,X] 1 ) 2 + (d 1 x 2 + A^ 
+ (d 1 x 3 - A\x 2 ) 2 + (d 2 x 2 + Alx 3 ) 2 + (9 2 x 3 + A\x 2 ) 2 

+ ^ 2 ) 2 + (* 3 ) 2 -!)} (31) 
For the case of A = \ this can be brought into the form: 

E = IJd^iiBlf + iBlf + iiA^XYf + HA^XY) 2 

+(Bl ± i((x 2 ) 2 + (x 3 ) 2 - l)) 2 + {{d lX 2 + A\x 3 ) t (d 2 x 3 - A\x 2 )) 2 
+((d 2 x 2 + A\x 3 ) ± (d lX 3 - A\x 2 )) 2 ± [A 1 ,A 2 ] 1 ((x 2 ) 2 + (x 3 ) 2 - 1)} 
±1 J rf 2 x( J B z + [A 1 ,A 2 ]) 1 (32) 

Still, there appears to be a problem, because there are two terms in this expression 
of the energy density, which are not squares and as we will show later, only one of 
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them is proportional to the winding number. The other term, [A ± , A 2 \ 1 ((x 2 ) 2 + (x 3 ) 2 — 
1), therefore appears to be problematic. This problem can fortunately be cured rather 
straightforwardly. Remember that we are already in a 'gauge' A t = 0. where furthermore 
the static fields are time independent. In this situation the residual gauge freedom of 
time independent gauge transformations may be used to put the term [Ai, A2} 1 equal to 
zero. In this gauge the energy density consists only of squares and a term proportional to 
the winding number. The minimum of the energy is now easily obtained by putting all 
squares in the energy density equal to zero. This then yields a set of first order equations of 
which the solutions are also solutions of the full field equations. The first order equations, 
including the gauge conditions, are: 



[Ai,A 2 ] 1 


= 


(33) 


Bl 


= 


(34) 


Bl 


= 


(35) 


[AuX] 1 


= 


(36) 


[A^X] 1 


= 


(37) 


{Bl±\{{x 2 ) 2 + {x 3 ) 2 -l) 


= 


(38) 


(d lX 2 + A\x 3 ) t (d 2 x 3 - A\x 2 ) 


= 


(39) 


(d 2 x 2 + A\x 3 ) ± (d lX 3 - A\x 2 ) 


= 


(40) 



The last three equations are identical to those that were obtained in the NO model. The 
energy of solutions to this set of equations are fully determined by the term J d 2 x(B z + 
[Ai,/^]) 1 which is proportional to the winding number, as we show next. 

The general expression for X in the presence of an Alice flux in the first isospin direction 
"along" Ti, becomes: 

X = x(r)e 2wix(9)Tl T 2 = x(r) cos(2tv X (9))T 2 + x(r) sm(2n X (9))T 3 = aT 2 + bT 3 (41) 

With x(9 + 27r) = x(0) + §■ For r — > 00 one has x(r — > 00) = 1, and the winding number 
can be extracted from the asymptotics by: 

— (f> d ln(a + % * b) = n = - (42) 
2vr J 2 

For r — > 00 one also has the spatial covariant derivatives DX = or: 

OX = [A, X] (43) 

or in components: 

[AX] 1 =0 

da = -A x b 
Ob = A 1 a 

From this one finds <91n(a + i * b) = i * A 1 , which means that: 

■ rr — »oo -1 rr — ^oo -i p 

— <b d\n(a + i*b) = — d> dll-A 1 = — d 2 x(B 1 z + [A 1 ,A 2 } 1 ) (44) 
2n J 2tx J 2-k J 



10 



Thus the rescaled energy of the solutions is equal to | . Note that the above expressions 
do not look gauge invariant because we are evaluating a gauge invariant expression in a 
particular gauge. 

One should, of course, check whether the first order equations actually do have any 
solutions of the type we are interested in. By inserting the ansatz used before and putting 
c(r) to zero, one arrives at the following set of coupled non-linear first order equations. 

rd r a{r) = -(1 — a(r))a(r) (45) 

-d r a(r) = l-a 2 (r) (46) 
r 

Now, these turn out to be a special case of the equations encountered before by De Vega 
and Schaposnik in their study of the NO model. They where obviously only interested 
in the case of integer winding number, whereas we are interested in the case of fractional 
winding number n — |. The corresponding numerical solution is given in figures and |3|. 

We have attained our goal of obtaining a set of first order equations, of which the solutions 
are also static minima of the energy (with no electric fields). As is well-known , first order 
equations play a deep role in gauge theories. Bogomol'nyi || explained, for the NO 
model, that solutions which come from the first order equations are also minima of the 
energy, which implied the neutral stability of such solutions. Later it was shown || that 
the occurrence so-called Bogomol'nyi equations is tightly connected to the existence of 
a super-symmetric extension of the theory. The explicit super symmetry extension of 
the NO model was given by |10| and in agreement with || showed that the first order 



equations indeed follow from an increase of super-symmetry. In our models we also found 
first order equations whose solutions are also solutions to the full set of second order field 
equations. We showed that the solutions are also minima of the energy. This obviously 
raises the question if these first order solutions can also be explained by an increase of 
super symmetry of a super-symmetric extension of our models. A superficial analysis 
suggests that this is not the case, basically because we can only recover the Bogomol'nyi 
argument within the context of a very restrictive ansatz. In this respect the situation is 



similar to that encountered in the study of regular Z^r monopoles JTTI] 



3 Conclusion 

In this paper we proposed two new models which both possess an Alice electrodynamics 
phase. For both models we constructed solutions corresponding to a topologically stable 
Alice flux. We found a way to project the theories on the Nielsen Olesen model and, in that 
way, obtained first order equations. Solutions to these first order equations corresponding 
to minima of the energy (without electric fields) were constructed numerically. 



We close with a brief remark concerning the zero modes of our solution. E.Weinberg [|12 
showed that in the NO model, for the critical value A = |, a flux with winding number n 
has 2n zero modes. These modes are interpreted as being the positions of the unit fluxes. 
At first sight this appears to give problems for the case of n = |, but carefully redoing 
section IV of the article mentioned, in particular using the fact that our fields are allowed 
to be double valued, one may show that the answer for n = | is that there are again two 
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zero modes, as one would expect. 

We thank M.M.H.Postma || for his contributions in the early stages of the project. 
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